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Abstract 

We compute the low energy mass spectrum of lattice QCD in the large Nq limit. Expanding 
around a gauge-invariant ground state, which spontaneously breaks the discrete chiral symmetry 
we derive an improved strong-coupling expansion and evaluate, for any value of Nc, the masses of 
the low-lying states in the unflavored meson spectrum. We then take the 't Hooft limit by rescaling 
g 2 Nc — > g 2 ; the 't Hooft limit is smooth and no arbitrary parameters are needed. We find, already 
at the fourth order of the strong coupling perturbation theory, a very good agreement between the 
results of our lattice computation and the known continuum values. 
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I. INTRODUCTION 

Since the seminal work of 't Hooft the large Nc limit (with g 2 Nc fixed) has played 
an increasingly important role in studying gauge theories in the continuum and on the 
lattice and, more recently, the duality between gauge and string theories jj, O] . Our aim 
is to use the 't Hooft limit to investigate some features of the meson spectrum of strongly 
coupled lattice QCD. Our results exhibit a good agreement with the observed meson masses, 
thus implying that in this context also the 't Hooft limit offers a very accurate method to 
analyze the QCD spectrum. 

A quite successful approach to studying gauge theories with confining spectra is the 
strong-coupling expansion In the strong coupling limit, confinement is explicit, the con- 
fining string is a stable object 7], and some other qualitative features of the spectrum are 
easily obtained. The formulation of the strong coupling expansion requires a gauge invari- 
ant ultraviolet cutoff, which is most conveniently implemented using lattice regularization. 
One of the most difficult problems of the strong-coupling approach to lattice gauge theory 
remains its extrapolation to the continuum limit, which usually occurs at weak coupling. 
In spite of this difficulty, there are strong-coupling computations that claim a high degree 
of success. A useful test of strong-coupling expansions in lattice gauge theories has been 
the study of lower dimensional models whose solution in the continuum is known even in 
the stron g co upling regime; for the Schwinger models [3] and the two-dimensional 't Hooft 



m „de. QU the m es„» S peet rUm a nd the ehha, eoadeasate have heen evaluated USlng the 
strong-coupling approximation and a remarkable agreement with the known exact results 



was found. Furthermore, it is by now well known [ll| that relevant features of strongly 
coupled lattice gauge theories have analogues in quantum spin systems: in many instances 
selecting the appropriate ground state for building the strong-coupling expansion of a gauge 
model turns out to be equivalent to finding the ground state of a generalized quantum 
antiferromagnet Q]. This equivalence turns out to provide pertinent nonperturbative in- 
formation about the structure of the gauge invariant states and, in many instances, greatly 
simplifies the computation of the chiral condensate. 

Many choices of strong-coupling theory produce identical continuum physics. As it is well 
known, any irrelevant operator may be added to the lattice Hamiltonian with an arbitrary 
parameter. Even if such ad hoc terms do not modify the continuum limit their coefficients 
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might appear in physical quantities such as the hadron masses Our study evidences 

that the 't Hooft limit does not require any arbitrary parameter and leads naturally to 
unambiguous and well defined results. 

We shall develop a method enabling us to extend to a general number of colors the 
ce.ebrated approach to the computation of the .neson spectrum in strongly coup.ed .att.ee 
QCD first introduced by Banks et al. in Ref. [13 . We shall use the Hamiltonian approach 
to lattice gauge theory using staggered fermions 1J|. The staggered formalism is known to 
yield good results in the strong-coupling evaluation of the hadron spectrum [l^; in contrast, 
other types of lattice fermions such as domain-wall or overlap fermions at strong coupling 
are expected to suffer both doubling and explicit breaking of chiral symmetry [15[. In the 
staggered fermion formalism, due to fermion doubling, the number of continuum fermions 
with Nf lattice fermions in D space-time dimensions is Nf2^ D ^ 2 ^~ 1 . If, as in this paper, one 
includes only one lattice fermion, the continuum limit yields the two lightest quarks (u,d). 
In this paper we shall determine the spectrum of the low-lying unflavored mesons. 

The computations are performed using x = 1/g 2 as the expansion parameter and the 
extrapolation to the continuum limit is carried by means of Pade approximants. As the 
ground state we take the one of the pertinent antiferromagnetic Ising model. The gauge 
invariant eigenstates of the unperturbed Hamiltonian are used to develop the perturbative 
expansion. The 't Hooft limit is then taken: the coupling constant is rescaled according to 
g 2 Nc — > g 2 and then Nc is sent to infinity. As shown in mesons for large N c are free, 
stable and non-interacting and their masses have smooth limits. Zweig's rule is exact at large 
Nc, mixing of mesons with glue states are suppressed and mesons for large Nc are pure 
qq states. The meson energies computed perturbatively in the strong-coupling expansion 
contain terms depending on the number of lattice links. In the thermodynamic limit these 
cancel against the ground state vacuum energy. Thus the results we shall present are finite, 
parameter independent and already in good agreement with the known physical values at 
the fourth order in the perturbative strong-coupling expansion. 

In Sec. II we review the well known Hamiltonian formulation of lattice QCD with stag- 
gered fermions and provide a classification of the symmetries of the theory. 

In Sec. Ill we consider the strong-coupling limit of lattice QCD and determine the chiral 
symmetry breaking ground state. We then evaluate its energy up to the fourth order in the 
strong-coupling expansion. 
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In Sec. IV we construct the operators creating mesons from the vacuum and then compute 
their energies up to the fourth order in the perturbative expansion. We show that these 
energies are finite and well defined in the large Nc limit. 

Section V is devoted to the extrapolation of the lattice results to the continuum theory. 
There we show that the lattice theory, properly extrapolated to the continuum via Pade 
approximants, already well reproduces the known experimental values for the ratios between 
meson masses at low orders in the strong-coupling expansion. 

Section VI is devoted to some concluding remarks, while the appendices illustrate some 
technical aspects needed to clarify the computations reported in the paper. 



II. LATTICE QCD IN THE HAMILTONIAN FORMULATION 

In the Hamiltonian formulation of lattice QCD with staggered fermions ^1 time is a 
continuous variable and space is discretized on a 3-dimensional cubic lattice with M sites, 
labeled by r = (x,y,z); with x,y and z integers. The conventions and the notation used in 
our paper are briefly reviewed in Appendix A. 

The lattice QCD Hamiltonian with one lattice flavor of massless quarks may be written 
as the sum of three contributions 



H = H e + H q + H m: (1) 



where 



H e = f a J2E a [?,nf (2) 

[r,h] 

H q = ^ E v(n)^ A (r + n)U AB [f, h]* B (f) + H.c. =H q + H\ (3) 

ZtiX r-* * i 
[r,n\ 

if m = — !— V \Tr(U\f,h}U\f + h,rh]U^[f + m,n]W[f,m\) + H.c.} (4) 

° [r,n,m\ 

are the electric field Hamiltonian, the interaction Hamiltonian between quarks and gauge 
fields and the magnetic Hamiltonian, respectively. The sums X)rf are extended to the N 
lattice links, whereas J2[?,h,m] is a sum over t ne plaquettes. n = x, y, z is the unit vector in 
the n direction and 

V (x) = rj(y) = (-1)*, rj(z) = {-If (5) 
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are the Dirac a matrices for staggered fermions 17 1. 

The gauge field U[f,n] is associated with the link [f,n] and it is a group element in 
the fundamental representation of SU(N C ). Two gauge fields occupying the same link are 
related by 

U[f, n] = W [f + n, — n\. (6) 
The electric field operator E a [r, n] is defined on a link and it obeys the algebra 



E a [r,h},E b [r',m] = if abc E c [f,h]5([f,h] - [r',m]) 



(7) 



E[r,n] = E a [r,h]T a , where T a , a = 1, ■■■,Nq — 1, are the generators of the Lie algebra of 
U(Nc)- They generate the left-action of the Lie algebra on U[f,n] 12. Il8] 



E a [f, n],U[P, rh) = -T a U[f, h)5([r, h) - [P, m\) 
E a [f,n},U^[r',rh] \ = U ] [r,h}T a 5({r,h] - [P,m]). 



(8) 
(9) 



The fermion fields \1/ are defined on the lattice sites and obey the anticommutation relations 

{* A (r), tf B (r')} = <W(r - P) (10) 
{V A (r),* B (P)} = 

In addition to gauge invariance,the lattice Hamiltonian is invariant under the following 
discrete symmetries 



1. Lattice translation by even integers, 



^(x, y, z) -> ^(x + 2l,y + 2p, z + 2q) 



(11) 



where I, p, q are integers. This symmetry operation amounts to a discrete translational 
invariance on the lattice. 



2. Lattice translation by a single link, 

*(r) *(r + x)(-l) y 
^(r) -> ^(r + y)(-l) z 
tf(r) -> #(r + z)(-l) x . 

In momentum space the last equation can be written as 



(12) 



q -> e ik *-y 5 r 3 q 



which in the continuum limit, where k z is infinitesimal, gives 



the other two transformations yield 



q -> 75^2? 
q -> 75 rig. 

The invariance under translation by a single link plays then the role of a discrete 
chiral invariance of the theory. This symmetry is broken by an explicit mass term in 
the Hamiltonian. As we shall show later, the chiral symmetry on the lattice is also 
spontaneously broken by the vacuum, and this should generate a nonvanishing chiral 
condensate < $ip >. 

3. Shift along a face diagonal, 

*(r) -> (-l) x+v V(r + x + z) (13) 
tf(r) -> (-l) w+z tf(r + £ + x) 
*(r) -> (-l) 2+a; ^(r + z + 

These transformations correspond to the discrete isospin rotations 

q -> r 2 g 
9 -»■ ^3? 
9 -»■ nq. 

4. Parity, 

*(r) -> *(-r). (14) 
This is the reflection through the origin. 

5. G-parity 

*(r) -> # f (-r). (15) 

This is just complex conjugation. 
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III. THE STRONG-COUPLING LIMIT AND THE VACUUM STATE ENERGY 



In this section we shall describe in some detail the method we propose for implementing 
the strong-coupling expansion for QCD with a generic number of colors Nc, by explicitly 
computing the vacuum state energy up to the fourth order in the strong coupling expansion. 

In the strong-coupling expansion the electric field Hamiltonian H e , Eq.(J2J), is the unper- 
turbed Hamiltonian while the interaction Hamiltonian between quarks and gauge fields H q , 
Eq.Q, and the magnetic Hamiltonian H m , Eq.Q), are treated as perturbations. Further- 
more, each term in H, Eq.Q, is gauge invariant since 

[^ a (f),^] = 0, i = e,q,m. (16) 

In Eq.(02|) 

n 

g\r) = £ E a [f,i] + & a (?)TZb9bW (17) 

i=—n 

are the generators of static gauge transformations and a — 1, ...,Nq — 1; these generators 
obey the Lie algebra 

[g a (f),g\r')] = if abc g c (r)5(r-f>). (18) 

The empty vacuum |0 > is a gauge-invariant ({? a (r)|0 >= 0) singlet of the electric field 
algebra (E a [f, i]\0 >= 0) which contains no fermions (^(r^O >= 0). 

Due to gauge invariance, |0 > must also be color singlet and should be charge neutral, i.e., 
it should obey the equation Y^fP{^)\^ >— 0, with p(r) given by the local fermion number 
operator 

p(f) = -[^(f), *(r)] = * f (r)*(r) - — . (19) 

At a given site one may form a color singlet in two ways: one can leave the site unoccupied or 
one may put on it iV<7 fermions with antisymmetrized singlet wave function. Fermi statistics 
allow at most one singlet per site and it is possible to distribute the M/2 singlets arbitrarily 
(M is the total number of lattice sites). The degeneracy of the lowest eigenstate of H e is 
then given by 

Ml 




The degeneracy may be removed by diagonalizing the first nontrivial order in perturba- 
tion theory. All matrix elements in the vector space of degenerate vacua of the first order 
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Hamiltonian vanish, since < 0|i7 g |0 >= (as / dl/U^B = E^l)- The first nonvanishing 
contribution comes from second order. At this order one should compute 

E® =< 0\H m \0 > + < 0\H n ° H q \0 >, (20) 

tin — ri„ 



where 



<, >= 



II / dU[r 



is the inner product in the full Hilbert space of the model. dU is the Haar measure on the 
gauge group manifold and (, ) the fermion Fock space inner product; n is the projection 
operator projecting onto states orthogonal to |0 >. In Eq. (|2()jl Ho is ineffective since the 
states created when H q acts on the vacuum are always orthogonal to |0 >. 

(2) 

Eq may be more conveniently computed by constructing an eigenstate of H e and using 
it to evaluate the function f(H e ) = E ^° H appearing in Eq.(|2U|l. Since the vacuum state 
|0 > is a singlet of the electric field algebra, one has 

E a [r,h]\0 >= (21) 

which, in turn, implies that 

# e |0>=0. (22) 

Using then Eq.(jHJ) and Eq. (J22J) and putting the commutator [H e , U[r,n]]\0 > in place of 
H e U[f,n]\0 >, one finds 

H e U[r,n}\0 >= ^C 2 {N c )U[r, n}\0 >, (23) 

where Ci = (Nq — l)/2Nc is the Casimir operator of SU(Nc)- U[f,n]\0 > is then an 
eigenstate of H e with eigenvalue g 2 C 2 (N c )/2a. Consequently, 

< 0\H q —±—H q \0 >= _-£L < 0|flJ^,|0 > . (24) 

After integration over the link variables U [19] [see also Eqs. (|C3J) . (|C4J) ] one sees immediately 
that the first term in Eq.()20|) is zero and that the second order correction to the vacuum 
energy is given by 



N c 

p(r + h) + — 



N c 

-P{r) + — 



|0 > . (25) 



S 



The diagonalized effective Hamiltonian is, up to an additive constant, given by 

H -=^brc^ w+AMif > 1 (26) 

a [r,n\ 

In deriving Eq. ()2fi|) one should also take into account the fact we that gauge invariant states 
such as |0 > should be charge neutral, i.e., X^p(r) = 0. As it is well known, Heft takes 
the form of the Hamiltonian of an antiferromagnet of spin Nc/2. It was shown in that 
even for domain-wall fermions the effective Hamiltonian is that of an antiferromagnet but 
the fermions are massive and doubled. 

Since, due to Eq. (fT9*j) . p(r) has only two possible eigenvalues p = ±N c /2, the lowest 
energy a link can have may occur only when one end has p = +Nc/2 and the other has 
p = —Nc/2. In the space of pure fermion states the true ground state must thus minimize 
H e Q and be fluxless. 

There are two ground states of H e g corresponding to those of the antiferromagnetic 
Ising model of spin Nc/2. In the large Nc limit the two ground states are not mixed at any 
finite order of perturbation theory. In fact, it would be necessary to apply H q at least Nc 
times in order to transform one ground state into the other, since H q acts as an hopping 
Hamiltonian destroying a quark on a site and creating it on a neighboring site. One may 
choose as the ground state the one in which p = +Nc/2 on even sites and p = —Nc/2 on 
odd sites, the other state being obtained by interchanging odd and even sites. With this 
choice of the ground state the sum over the lattice links [f, n] in Eq. (|25|) may be easily done 
and, for the ground state energy at the second order in the strong-coupling perturbative 
expansion, one gets 

^ - ~^k N - (27 > 

As evidenced in Eq.()12jh the chiral symmetry on the lattice is given by the translation by 
a single link and takes even sites into odd sites. Having chosen one of the two ground 
states described above, chiral symmetry is spontaneously broken in the large Nc limit. 
Thus, in the perturbative expansion one has to consider only diagonal matrix elements and, 
consequently, perturbation theory for nondegenerate states. 

The ground state energy up to the fourth order in the strong coupling expansion, is given 

by 

E = E^ + eP+E^\ (28) 
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where 



with 



-^0 4) = [ E O i ']l + II + [Eq 1 '}!!! 



,(4) 



(4)i 



[4% <0\H, 

n(4) 



n, 



9 E^ 



-H„ 



n n 





< 0\Hn 



H e q E^ 



-H, 



n f 



9 77(0) 



n r 



H q \Q >< 0\H q 



n 



H q \o> 



(4 0) - H e 



-HJ0 > 



-H m \0 > . 



(29) 

(30) 
(31) 
(32) 



Equation (J3~T|) has the same form of a second order contribution and may be evaluated 
following the same steps used to arrive at Eq.()27|): the only difference being that, in Eq.()3H). 
the energy denominator also appears squared so that, using Eq. (|2^j) . one gets 



< 0\H n 



1 



H n \0 >-- 



(E — H e j y 1^2 

After integration over the link variables, one finds that 

n fn„. _ lfr n 



< 0\HlH a \Q > ■ 



g4C 2 - «l 



(33) 



< 0\H. 



&o — ti e 



Hg\0 X 0\Hg~ 



(4 0) - H e 



-HJ0 >-- 



g«aClN% 



< 0| J2 n (?+ n){n{f) - N c )\0 > 



N 2 C N 2 
2g 6 aC$ : 



where 



n r 



p(r) + 



N c 



(34) 



(35) 



yields the number of fermions at the site r. 

One may now turn to the first term contributing to Eq 4 \ Eq. (j30|) . Using Eq. (j23)l . one 

gets 

n ~ n ~ n 



< o\h, 

4a 2 



9=7(0) 
^0 



# g |0 >= 



(36) 



4 < OliLif! 



n 



q i E (o) 







fl^lO > +2 < 0\H o H, 



Tic 







> 



In deriving Eq. (|3fi|) from Eq. (j30j) . one should observe that the most external Ilo is ineffective; 
the only task to accomplish is then to evaluate the energy denominator in the middle. For 
this purpose it is most convenient to rewrite H q H^\0 > and H^H^\0 > as linear combinations 
of eigenvectors of H e , and in order to do this one should first consider the action of H e on 
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UW at different links. Using the commutators (JHJ), @ and the fact that H e \0 >= one 
finds 

H e UAB[r,n}Ul D [?,m]\0> = ^ (C 2 (N c )U AB [r,n}Ul D [P,m] 

+ 7^T U ab [r, n}U ] CD [f , n}8( [f, n] - [P, rh] ) 



l5 AD 5 BC 5{[r,n] - [r',fh])) |0 > • (37) 



Keeping into account the action of n it is easy to see that U AB [r,n]Ul< D [r\rh}\0 > is an 
eigenstate of f(H e ) = U /(Eq ^ — H e ) with eigenvalue 



f(H e )U AB [f,n]U^ CD f',rh}\0 >= -—U AB [r, n]tf CD [r', m}\0 > . (38) 

A similar procedure shows that when H e acts on two U's at different links, one gets 

H e U AB [r,n]Ul D [r>,m]\0> = 9 - (C 2 (N c )U AB [f,h}Uh D [?,rn} 

~ ^^U AB [rMUh D [rMK[rM-[r'M))\^> ■ (39) 

In Appendix B we show that the combination (u AB [r, n\U}j D [r\ rh] — U AB [f, n]Uc D [r, n\) |0 > 
is indeed an eigenstate of H e with eigenvalue g 2 C 2 (Nc)/a. Using this result and the fact 
that U AB [f. , n] Uq D [f. , ft] | > is an eigenstate of H e one can easily show that 



f(H e )U AB [f,h]U' OD [r',m]\0 >= -— (U AB [r, ft]Ul D [r>, fh] 

+ A^2 f/ ^ [f '^^ D[f '^ 5([f '^ " [ ^™ ]) ) 10 > ■ (40) 
Taking into account Eq. (jHHj) and Eq.pH|). Eq. (j3T)j) becomes 

< ^V ^ H Hq F^° H Hq F^ H Hgl ° >= 
'"" (-4 < 0\H q HlH q Hl\0 > -2 < 0|ff,#,ffjffj|0 > 



1 < 0|iy gJ Ff 9 ^^(r)C/i B [f,n]* B (r + n) 



2a 2 (iV c -2) 

#^tf<yr,n]^(r + ra)|0>). (41) 



Integrating now over the link variables (see Eqs. ()C4|) . (jC5j) ). one gets 

< ®\H<1 — (W^ H q — 77y— H q — 77y-^ Hq\0 > = 

En — H P Ek°>-H P En — H P 
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1 



'g«aCi 



1 



< ^1 I ~~ Z n(f+ h)(n(r) — Nc)(n(r + h + rh) — Nc) 
+ E n(f + h)(n{f) — N c )n{f — m) + E n(f + n)(n(f) — Nc)n(f+ rh) 

[r,h]^[r— rh,rh] [r,n\^[r;m\ 

E n(r*+ h){n{f) — Nc)(n(f + n — m) — iVc) 

[r, n] ^ [f -I- h— rh , rh] 

+ E n(r + n){n{r) — Nc)n{r' + m)(n(r') — iVc) 

[r,n]^[r',m] 



+ E ra(f + h)(n(r) - N c )n(r')(n(r' + m) - N c )j |0 > 

< °l E (-n(r + h)(n(r) - N c ) + n(f + h) 2 (n(r) - N c ) 



[r,h\^[r',rh\ 
1 



i\r c (AT c -2) 

-n(f + h)(n(f) - iV c ) 2 + n(f + h) 2 (n(r) - N c ) 2 ) |0 > 



2^6aC 2 3 



2 A72 , iV c (10iVc - 21) 



N, 



N 



c 



(42) 



In deriving Eq. (|42|) one should use Eq. ()C5|) when two C/'s and two L^'s are on the same link. 

Equation (J32|) in is the magnetic contribution to the vacuum energy. Since H m 
contains already a factor 1/g 2 , it contributes only at the fourth order of the strong coupling 
expansion. It can be written as 

n 



< 0\H, 



m E (o) 



H, 



H m \0 >-- 



4<0 I E TTT" u ab [r , n] U BC [f + n, rh] U CD [f + rh , ft] A [r. , rh] , Q) n ° 



^ 2# 2 a 



&o — n e 



E 1T2-U EF [?, l]U FG [r> + 1, k]U GH [? + k, l]U j HE {P, k]\0 > 
,tt~ lira 

[r',l,k] ^ 

Using Eq.(jHJ) and Eq.®, one finds 



(43) 



H e U AB [r. ; ft] U BC [r + n, rh] U CD [f + rh, ft] U^ DA [f. , rh] | >= 

2 ^C 2 C^b [f ; n] £/ BC [r + ft, rh] U CD [r + m, ft] U ] DA [r , m] . (44) 

Thus, a plaquette acting on the vacuum |0 >, Tr[UUWW]\0 >, is an eigenstate of H e with 
eigenvalue 2g 2 C 2 /a. Equation (J4H)) then becomes 

n 



<0\H n 



- H, 



H m \0 >- 



E ^[r', f]^ G [P + f, £^[r' + fe, f]^L[r', fc]|0 > . 
■?,i,k] 



(45) 
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The integration over the link variables requires that the two plaquettes are coincident, oth- 
erwise the integral would vanish. Keeping into account that 2N is the number of oriented 
plaquettes on the lattice, one gets 

< 0\H m ^^ H m \0 >= --^r. (46) 

Collecting all the terms ()34j) . (J42j) and (J46|) the fourth order correction to the vacuum 
energy becomes 



(4) _ iv iv Cl xuivc - ±±) iv _ {4 ) (4) . 

L ° ~ 2^Cl N^2 7^2 = q ° m °' { } 

where E$ is the contribution due to the quark Hamiltonian and E^l the one due to the 
magnetic Hamiltonian. Note that is proportional to N and thus the vacuum energy 



is an extensive variable. Moreover, the N 2 dependence of Eq. (134)1 is precisely canceled by 
Eq.(|42j) . As a check of our computation one may set Nc = 3 in Eq.([47|) and compare the 
result with those obtained in d3] using a completely different approach. The agreement is 
exact when the coefficient of the irrelevant operator introduced in is set to zero j^ . 



IV. THE LOW-MASS MESON STATES 



In the strong coupling expansion the lowest-lying states in the meson spectrum are those 
consisting of a quark and an antiquark at opposite ends of a single link. If the quark is at 
(r + h) and the antiquark at r a basis for such states is given by 

\r,n>=&(r + n)U[r,n]iff(r)\0>. (48) 

For a given meson, the wave function may be determined through the following steps. One 
may first take the quark bilinear in the continuum with the desired transformation properties 
and appropriate continuum quantum numbers and then write it in point-separated lattice 
form using the discrete symmetries of the theory. Only after fixing the pertinent lattice 
quantum numbers one may apply the bilinear to the vacuum. 

In the following we shall be interested only in the low- lying unflavored mesons: ttq, p, u, 
bi, en-, f 2 Jo- In the continuum theory the wave functions for these mesons are given by 

|tt >~ ?"*75^r 3 ^|0 > (49) 
13 



\uj >~ ^W^lO > 



1 



\p >~ Wa x -r 3 V\0 > 
\b x >~ i^j 5 d z r 3 ^\0 > 
\ ai >~ i^(a x d y - a y d x )T 3 ^\0 > 
| /a >~ itf^a^ + a s 5 x - 2a y d y )V\0 > 

\f >~ y^|o > . 



(50) 
(51) 
(52) 
(53) 
(54) 
(55) 



The choice done is as in and it is based on the quantum numbers labeling the mesonic 
states. For example, for 7r one needs a pseudoscalar with nontrivial isospin; thus Wj^T^ty 
is a pertinent choice of the wave function. The choice for the components of the vector 
mesons (u, p, b%, ax) or of the spin-2 meson (/a) is made by observing that these are the 
only components of these mesons that on the lattice have the standard form f)48j) . 

The lattice form of these operators may be obtained by applying the staggered fermion 
formalism to derive operators with appropriate lattice quantum numbers. The lattice wave 
functions at zero momentum are then given by 



Fo > 
\uj > 

\P > 
\bi> 
\ai > 

1/2 > 

l/o> 



VNcM 
i 

\fNcM 
1 

1 

i 



Y l (-l) x & A (T f + z)U AB [r,z}V B (r) - h.c. 

r 

£(-l)^(f + z)U AB [r, z}V B (r) - h.c. 

r 

J2(-i) y ^A(r + y)U A B[r,y}^ B (r)+h.c. 

f 

£(-ir^(f + z)Uab[t, z}V B (r) + h.c. 

r 

Y^{-l) y ^\{r + v)UAB[r,m B {f) 



|0 > 
|0 > 
|0 > 
10 > 



+ {-l) x+y+z ¥ A {r + x)U AB [r,x]y B (r)} - h.c 



10 > 



Y^[(-l) y ^ A (f+z)UAB[r,z]^ B ( 



V2N C M 

+ (-l) z ¥ A (f + x)U AB [f, x]V B (r)] + h.c. 



(56) 
(57) 
(58) 
(59) 
(60) 

(61) 



10 > 



y^y(n)W A (r + n)U AB [f,n]^ B (f) + h.c. 



10 > 



y/3N c M 

The normalizations are fixed in the standard way by integrating over the link variables. 
All the mesons are degenerate at the lowest order and their energy is given by 

E$ =<M\H e \M > 



(62) 



9 2 C 2 



2a 



(63) 
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as it can be easily seen using Eq. (j23J) and integrating over the link variables using Eq. (jC4|) . 
Since a static limit (fixed and large a) is used -and thus the states cannot propagate in this 
approximation- all the single-link mesons have the same mass regardless of the character 
(e.g., s-wave or p-wave) of their continuum wave functions. As we shall see, the fourth order 
computation will cure this unphysical effect of the static lattice approximation. 

The meson energy has been computed up to the fourth order in the perturbative expansion 
using a method similar to the one used in the evaluation of the vacuum energy. If one 
evaluates the meson energy up to the fourth order, one gets 

p _ 9 2 C 2 p (2) p (4) 



where 



and 



E$ =< M\H q , Q nM H q \M > (64) 



P» TT 



<= <M\H q ^ F H q ^ ¥ H q ^ ¥ H,\M> (65) 

- < ^ft^r^ ><: M ^ ( ^ M „,M M > ^ 

+ <M\H m " M H m \M>. (67) 

Ha/ is the projection operator projecting onto states orthogonal to those states that have 
the unperturbed energy of a generic single link meson (|6Hjl . \M. > is the meson operator. 



A. Second order 

The projection operator Hm does not affect the second order in the strong coupling 
expansion since the states created by Ai and H q acting on the vacuum |0 > are orthogonal 
to the meson state. Thus, the matrix elements to be computed are 

E$ =<(M + M*)\(H q + Hi) 1 (H q + H\)\{M + M*) >, (68) 

where \M. > is one of the meson operator (|56|) - (j62|l and M. = M + . In Eq. ()68)1 the terms 
containing a different number of U and W vanish when integrated over the link variables. 
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The remaining terms may be conveniently grouped as 

E$ = 2 [< M\H q K M {H e )H\\M^ > + < M\H\K M {H e )H q \M^ > + < M\H\k M {H e )H\\M > 

(69) 

where 



E {0) - H ' 



All the one-link meson states ()56jl -(|62 jl consist of linear combinations of directed links on 
the lattice with appropriate phases, which are responsible for the differences in the matrix 
elements of the various mesons (|56| ) -([62 p . 

To elucidate the method used to compute these matrix elements, it is sufficient to consider 
a generic meson of the form 

^S(^ A (r + n)U AB [r,n}V B (r)\0 >; (70) 

f 

In Eq. ()70p S(r) is one of the phases of the meson operators (|56| ) -([62 p . In order to compute 

(I I) i =< M\H q K M {H e )H\\M^ > (71) 
(J/) 2 =< M\Hl~K M {H e )H q \M^ > (72) 
(//) 3 =< M\H\K M {H e )Hl\M > (73) 

one should construct suitable eigenstates of the unperturbed Hamiltonian H e in order to 
evaluate the function of H e , A M (H e ). Using the eigenvectors found in the previous section, 
it is quite easy to obtain 

A M (H e )U AB [f,n]Ul D [?,m]\0 >= 
2(1 [-U\ B [r, h]Ui D [P, m] - j^^uUr, ft}tf CD [r, h]j |0 > (74) 



and 



9 2 C 2 



A M (H e )U AB [r,h]Ul D [P,m]\0 >= 

-^r ~Uab[t, h]U}, D [?, rh) + 2 - U AB [r, ft)U^ D [r, ft]6 ([f, ft] - [r',m]) 
g L>2 \ + 1 

+T^t8 ad 5 bc 6 ([r,n] - [P,m])) |0 > . (75) 



N 2 C + 1 

Plugging Eq.(J7U) and Eq.(J75J) in Eqs.JZU-JZSIl, one gets 



= - — - < M\H q H\\M^ > 

9 L/ 2 
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1 



< M\H q -J2v(nm(f f )U\ B [f,h}'9 B (f+ 



n 



g*C 2 (N c -3) 1 *2 lr(i] 

^S*(^U^U j CD [r,h]^f D (f+n)\0 > . (76) 



(H)2 = - < M|fl}fT fl |Jlft > 

1 < M l^jJ E *7(™M(r + n)C/ AB [r,n]* B (r) 



+ g 2 C 2 (N* + l) N1 "'"«2 



t 



E^(^c(^cd[^]* d (? +n)|0 > 



4iV c 



02^0 + L ) 2 [r-,n] 

]T 5*(r)^ ; (r)1' D (r + h)5 AD 6 BC \0 > . (77) 



(^)s = - T^T < ^IfljAjl^ > 



^,S(r)^(f+n)f/ CD [f,n]^ D (f)|0 > 

r 

+ <M[gtl^ry(n)^(r> B (r + n) 

J2 S(f)^ c (r + n)y D (r)8 AD 5 B c\0 > ■ (78) 



One now needs the results of Eq. (jC4|) and Eq. ()C5|) to integrate over the link variable U and 
to choose the particular S(r) selecting a given meson. 

As an example, consider the case of the p meson. In terms of the fermion number operator 
n(r) one has 

1 



l(//)ll ' = -2?^ <01 



} j n(r + h)(n(r) — Nc)(n(f + m) 

f,[r,y] 7^[r,m] 

~~ E n(r + y) (n(f^ — Nq) (n(r + y — rh) 
+ n(r + y)(n(f)-N c )n(R + m)(n(R)-N c ) 

r,[r,y]^[R,m] 

+ E {-l) x+y+X+Y n{r + y)(n(r) - N c )n{R + y)(n(R) - N c ) 

r,R,{r,y]^[R,y] 

IN 

+ Nc E {-n(f+y)(n(r) - N c ) - n(f+ y){n{f) - N c f 
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+n(f + y)\n{r) - N c ) + n(r + y) 2 {n{r) - N c ) 2 )] |0 
2g 2 aC 2 I Nc — 3 4 J 



> 



(79) 



n(r + y)(n(r) — N c ){n{f + y + m) — iV ( 



-VScyvT °l 

+ X] n(r + y)(n(r) — A r c)n(r — m) 

f,[r,y]y^ [r— m,m] 

+ 51 n(f+y)(n(f^-N c )n(R)(n(R + m)-N c ) 

r,[r,y]y^[R,m] 

+ E (-l) wX+y n(f+y)(n(f) -N c )n(R + y)(n(R) - N c ) 

r,R,[r,y]?[R,y] 

( N 2 

+ E ^2 y) 2 («(r) - Nc) 2 - n(f+ y)(n(r) - N c )\ 

- jV f^ 1 [-^(r + y) 2 (n(f) - jV c ) + n(f + y)(n(r) - Nc) 2 ]j 
N 2 ( 

[EHr + y)\n{r) - N c f - n(r + y)(n(r) - 7V c) ] 



N 2 c + l 



-j^rj Y\- n (? + £) 2 W*0 - + + £) Wr) - N c) 2 }) 

+ E (-ir +?/+x+y n(f +y)(n(f) -N c )n(R + y)(n(R) - N c ) 

r,R,{r,y}^{R,y] 

?N 2 1 
-^fYE(-ir +?/+X+ Mr + y)(n(f) - N c )n(R + y)(n(R) - N c ) f |0 > 



\R 



-J—U-^N) 

2g 2 aC 2 [ 6 4 7V ' ' 



(80) 



1 



<0| 



j3Jp 2(? 2 aC 2 iVg 
- E^+ - iV c )(n(r + 2y) - 7V C ) 



s ^2 / n{r + y){n{f) — Nc)n(f — y) 

1 



|0>= 



2g 2 aC 2 



Adding up the three terms, one finds 

= 1 

p 9 2 aC 2 



' 2_(Nc-l)_Nc N 
N c -3 2 



(81) 



(82) 



A similar procedure yields the following results for the second order correction to the 
meson energies 



= 

*° g 2 aC 2 



' 2jNc-l 1 _Nc N 
N c -3 2 



(83) 
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1 

g 2 aC 2 
1 

g 2 aC 2 
ai 9 2 aC 2 



E, 



(2) 

h 



E 



(2) 
h 



E 



(2) 
fo 



g 2 aC 2 
1 

g 2 aC 2 



5 - 
7- 
9 - 
9 - 
11 - 



2(A C 


-1) 


N c 

-N 

2 


Nc 


-3 


2(N C 


-1) 


N c 

—N 

2 


N c 


- 3 


2(N C 


-1) 


N c 

—N 

2 


N c 


-3 


2(N C 


-1) 


N c 

-N 

2 


N c 


-3 



2(N, 



c 



N 



c 



N 



(84) 
(85) 
(86) 
(87) 
(88) 



N c -3 2 

where the A-dependent terms cancel against the vacuum energy, as should be since the 
masses are intensive quantities. 



After rescaling the coupling constant according to the 't Hooft prescription (large Nc 
with g 2 N c fixed) 

g 2 N c - g 2 

one finds that the large Nc limit makes Eqs. ([82 )1 - ([88 )1 finite. 

Our results have been checked by deriving the second order meson masses for Nq generic 
with the graphical procedure of Ref.|3|. In this case also there is complete agreement. 



B. A comment on irrelevant operators 

Since the strong-coupling limit of QCD is not universal, adding an irrelevant operator to 
the Hamiltonian, leads to the same physical predictions in the continuum limit. This allows 
one to introduce arbitrary parameters, the coefficients of these irrelevant operators, which 
are then fixed by fitting the experimental data. Our analysis shows that in the 't Hooft limit 
(Nc — > oo, g 2 Nc fixed) the meson masses can be made independent of arbitrary parame- 
ters and that results in agreement with experiments can be obtained without introducing 
irrelevant operators. 

In the celebrated computation of the hadron spectrum by Banks et al. Q| the lattice 
Hamiltonian was indeed modified by the addition of an irrelevant operator given by 

W = A]T (p(r)p(f+ n) + ^ j (89) 
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where A is a dimensionless irrelevant parameter. The new term was chosen according to 
three demands. First, it must remove the degeneracies at zeroth order so that nondegenerate 
perturbation theory can be used. Second, it must preserve the symmetries of the original 
Hamiltonian. The vacuum state of the modified theory must again break the chiral symmetry 
spontaneously. Third, the added term should have no effect on the continuum limit of the 
lattice theory, so it should be an irrelevant operator. W is a four-fermion operator and, when 
it is written in terms of the continuum variables with the conventional units, it depends on 
Aq 2 a 2 so that in the continuum limit (a — > 0) it vanishes faster than a 2 . 

n 

The irrelevant operator ()89j1 was introduced by Banks et al. mainly because in this 
way the meson masses are well defined even for Nq = 3 at the second order in the strong 
coupling expansion. The corrections ()82)) - (j88)) are in fact divergent for Nc = 3. This problem 
was avoided by introducing W in the unperturbed Hamiltonian. The meson masses then 
depend on the arbitrary parameter A. For Nc = 3 the diverging contribution, evidenced in 
Eqs. (182)1 - (}88|) . comes from a term in which the meson operators and the quark Hamiltonian 
are on the same link. In 13] this term varies as A -1 and if A is set to zero it yields a 
diverging contribution. It is not difficult to verify that, if the irrelevant operator is not 
introduced, up to the fourth order in the strong coupling expansion, there are divergences 
of the type 1/(N C - 1), 1/(N C - 2), 1/(N C - 3) and 1/(N C - 4); at the next order there 
are divergences up to \/{Nc — 5); and so on. In the infinite Nq limit these divergences are 
avoided. Our analysis thus shows that in the large Nc approach to lattice QCD there is no 
need for an irrelevant operator; in fact, with the 't Hooft prescription, the limit Nc — > oo 
yields series expansions for the meson masses which are free of divergences and thus well 
defined. 

The constant A in Eq. was fixed in jl^ by requiring that the qq state is less massive 
than a nucleon-antinucleon state in the static limit. As we will show below, in the 't 
Hooft limit however, the baryon masses are zero at zeroth order in the strong coupling 
expansion and acquire a mass proportional to Nc only at second order. Thus, baryons may 
be consistently regarded as QCD solitons and the unperturbed mass of a bound state 
nn in the large Nq limit vanishes at the lowest order in the perturbative expansion. A qq 
state is then not degenerate with a nn state and for this reason also there is no need of an 
irrelevant operator. 

Nucleon masses should be determined using an antisymmetric operator creating Nc 
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quarks at the same lattice site when acting on the vacuum state. The normalized nucleon 
state is 

\n >= ^j\/^E e A ^~ A "c tf M (f)tf M (f) . . . *i Wo (TO|0 > (90) 

At zeroth order in the strong-coupling expansion, the baryon is massless, since the creation 
operator \n > does not contain any color flux and thus H e \n >= 0. This is in agreement 
with the requirement that the nucleon mass, being the mass of a soliton, should vary as 
the inverse of the coupling constant. At the second order in the strong coupling expansion 
baryons already acquire mass given by 



E {2) =< n\H n 



n« ^ l 



N C 

-^N + N C 



m n — J^n ^0 — 2 r< 
g 2 aC 2 



N, 



c 



(91) 



which, after rescaling the coupling constant according to the 't Hooft prescription g 2 N, 



c 



g 2 , varies as Nc in the large Nc limit. 1/Nc is the "coupling constant"; thus the baryon 



mass again varies as the inverse of the coupling constant as a soliton mass should do 



1( 



Fourth order 



The fourth order corrections to the meson energies are given by the matrix elements 
(JE5J), (JEEJ) and (|57|). where the projection operator U M eliminates the states proportional to 
\M >. Again one can construct eigenstates of the unperturbed Hamiltonian to evaluate the 
function of H e in Eqs. ([fi5j) . (IFIHJl and <|B7 )l . 

First consider Eq. ()o5j1 . The non vanishing terms in Eq. ()o5j1 can be grouped as follows 

n M Hm Hm 



< M\H, 



E M - He 



P (0) TT Hq P (0) TT 



■H q \M > 



< | HlK M H\k M H q K M H q \ M > + < M\HlA M H\A M H q A M H q \M^ > 
+ < M\H q A M HlA M H q A M Hl\M^ > + < M\H\A M H q A M H\A M H q \Mi > 
+4 [< M\H q A M HlA M HlA M H q \M^ > + < M | H q Am H^Am H^AmH^ \ M > 
+ < M\HlA M H q A M H\A M Hl\M > 



(92) 



where 



E 



n 

w 

M 



M 
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To compute these matrix elements, first note that the two external projection operators do 
not affect the calculations. The projection operator in the middle, instead, does not allow 
patterns in which there are fermion operators creating and destroying the same quark at the 
same lattice site. For example, a term of the form 



^ A E ^A^)U AB [r,n}^ B {r + n)tf c (r + h)U CD [r,h}m D {r)\M 



'"e$-hM 



> 



is eliminated by the projection operator since it gives rise to a state with the same energy 
of a single link meson. 

In order to illustrate the method used in the computation of the meson energy, one should 
focus the attention on the first two terms in Eq.(|92|) and evaluate them for the generic meson 
([70)1 . Note that in the first of the two terms the projection operators are irrelevant. Let us 
define 

(A) =< M^\HlK M HlK M H q K M H q \M > (93) 

Using Eq.((7U), (J9*3~j) may be rewritten as 
4a 2 



(A) 



< M^\H\HlK M H q H q \M > 



N, 



c 



<M^\H\H\K M H t 



^(^^(f+njZ/ABt^nl^B^^^^r + nJt/ojjt^nl^ij^lO > 



<0| ]T S^f^rO^L^n^^ 



(N c - 3) 2 



rj(m)^ E (r' + m)U EF [r', m]^ F (f')S(r')^ G (r' + m)U GH {? ,m]^ H (r')\0 > 

r',[r' ,rh] 



(94) 



Making use of Eq.tJHJ), Eq.tjHJ), Eq. (J22j) and after constructing a suitable eigenstate of H e , one 
finds 



Am (# e ) Uab [f, n] U CD f', m] U EF [r 7 ', I] | > 
U ab [r, n] U CD [?, m] U EF [r 7 ', I] + 



g 2 C 2 • ' J ' N c -2 

+UAB[r,n]U C D[r',m]U EF [f, h] + U ab [r, ft] U C d [r' , m] U EF [r' , m 
6 



(UAB^npCD^flpEF^'.l] 



+ 



-U AB [f, n\U CD [r, n]U EF [f } n\ 



(N c -A)(N c -2) 
Taking into account Eq. (j93j) . Eq.fjHU) becomes 



> . 



(95) 
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An 3 

' {- < M^\H\H\H q H q \M > 



g 6 cl 

1 



N c -2 



< Aft|#}iftJL J2 +n)U GH [r,n\* H {7)*\{r + n)U IL [f,n]^ L (r)\M > 

[r,n] 

+ <M^\HlHl^ a ^r 1 {h)^{r + h)U GH ^ 

[r,n] f 

+ < M^\H\H\H q ±- Y,v{n)*W + n)UiL[r,n]^ L (r}J2S(^M(r + n)U MN [f,h]^ N (r)\0 > 



q q " H 2a^ 

[r,n 

6 



< M\H\Hl^- 2 ^Ur + m G H[r,n]^ H (^(f+n)U IL [f,n}^ L ( 



(Nc-4)(N c -2) 1 q q Aa\ 

J2S(mUr + n)U M N[r,h]^ N (r)\0 >\ 
f ) 

J2S(mli(r + n)U MN [r,h]^ N (r)\0 > 

] —— < M^\H\Hl £ *{,(r + h)U GH [f,h}^ G (^¥j(r + n)U IL [r,h}^ L (^ 



q <? 

\r,h 



a(N c - 4) 
J2S(mUr + n)UMN[r,n}^ N (r)\0 > 
4a(N c - 1) 



n 



g e C$(N c -3) 2 (N G -2) 

Rn] [r',rh] 

£ S(P)^ M (? + m)U MN [r', m]* N (?)\0 > 
P 

1 ___ t _ rrt 



< 0|$>*(fOn(0£4 J [^]Mr > + 



a(iV c - 4) 

WM [r',m] 

V] ?](m)\I / }(r / + rh)UiL 

[r',m] r' \ 

The matrix elements in Eq. ()96|) are evaluated by means of the integrals over the gauge group 
elements given in Appendix C. 

Group integration yields a nonvanishing result only if each link exhibits a combination 
of matrices from which a color singlet may be formed. Then, for each matrix element one 
should compute all the possible integrals obtained by putting three UW pairs on a different 
link [and then using Eq. ()C4j) ] or two of them on the same link [Eqs. (|C4J) and (|C5|) ]. or all 
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the three on the same link [Eq. (|C6j) ]. For p this term gives 



1 \ {17 

4 



4 



Xr | y + 4(iV c - l)(10JV o - 49) 5N C (N C - 1) 



(7V c -2)(AT c -4) 



2(N C - 2) 



iV 



iVc-3 
4(iV c - 



N C (N C -_1Y N A{N C - l) 2 (4iV c - 19) 



c 



(JV -2)(JV -4) 
2(iV%-8iV c + 22) " 
2(iY c -2)" {N C -2){N C -A) 



N, 



c 



-N 



(iVc-3) 2 
Now, one defines 

(5) =< M\H\k M H\k M H q k M H q \M^ > . 

Combining Eq.(|75j) and the action of Um, Eq.(|9*Hj) becomes 
4a 2 



(97) 



(98) 



(B) = ^ [< M\H\H\K M H q H q \M^ 



> 



16a 4 



— < M| X] ^(^^l B [r,n]^ s (f+n)^(fO^ D [f,n]* D (f+ri) 



A M «4 (r + n) C/ BF [r , n] * F (f) *k(r + h) U GH [r , n]^ H (f ) | M f > 



(99) 



Using again Eq.flBJ), Eq.® and Eq.(|22J), one finds 



A M (H e )U AB [r, h}U CD [r', m]C/^ F [r 7 ', Z]|0 >= 



2 2 Q 



■^Afltr, n]f/ CD [^,m]^ F [r 7 ', Z]|0 > 



(100) 



and 



A M {H e )U AB [r, n]U CD [r, n]U EF [r', 

- 9 \N c l a i)lN c -2) UAB ^ h]UcD[ ^ 
Taking into account Eq. (jl00|) and Eq. (|101|) . one has 

4a 3 



> 



(101) 



(B) 



< M\HlHlU M H q H q \M^ > 



+ 16^iV c -2) < M| ? *l(0^i B [^"]*B(r + n)^4D^^]*D(r + n) 



(102) 



U M ¥ E (f+ h)U EF [r, n]* F (r)V G (r + n)U GH [r, h}^ H (f)\M^ > 

The two different combinations of link variables allowed for the first term in Eq. (jl02|) are 
those where the gauge fields of the meson operators M and are defined on the same link. 
Then, the two possibilities may be represented by 



U AB [f, h]U CD [P, rh]U EF [?, l]U GH [?, l]U IL [?, m]Ul N [r, h] 
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U AB [f, n\U OD [r>, m\Ul F [?i l]U GH f', m]U IL [P', 1}U_ 



where the link variables in this expression have the same ordering they have in Eq. (jl02j) . 
Since the gauge variables of the quark Hamiltonian are on the same link, for the second 
term in Eq. (jl02|) there is only one possibility, i.e., the gauge fields of the meson operators 
must be defined on the same link, 



U A b [r. ; n] U CD [r ', rh] U\. F [r ', rh] U GH [r ', rh] U IL [r ', rh] U ] MN [r,n\. 



After the integration over the gauge fields one gets for p 
1 



4g 6 aC 3 2 



12(N C -1) 2 f Nl Nc(Nc-l) 2 

-72 - 12N C + ——^ + N 8N C + — ^— ^ f- 

6 N c -2 V 2 2 ( N c ~ 2 



N 2 

^N 2 
4 

(103) 



Analogously, the other terms of Eq. ()92j1 yield the following contributions for p: 

1 



< M\H q A M HlA M H q A M HW >-. 



N, 



C A7-2 



N z - 



-37- \2N C + N \AN C + " c 



2 



N, 



(N c -2)(N c -3yK U ^- 1) + 1 f^- 1)N 



< M\HlA M H q A M H\A M H q \M^ >-- 



-40 - 12iV c + 



+N 4N C + — ^j— r- —N 

V 2 2(N c -2) J 8 



< M\H q K M HlK M Hlk M H q \M^ > 



N C (N C - l) 3 \ 



-34 - 12iV c + 



( N 2 

+N AN r H - 

\ ° 2 2(N C - 3)(N C - 2) 



< M\H q k M H\k M H\k M Hl\M >-- 

< M\H q A M HlA M HlA M Hl\M >-- 



1 



N 2 
^N 2 



16 + 4N C - 



12(N C - l) 2 
N c -2 

12(iVc-l) 3 
(Nc-S)(N c -2) 

4(iV c - l) 3 



(N c -3)(N c -2) 



-N C N 



4(N C - l) 2 
16 + AN C - -\r - N cN 



N, 



c 



4g*aC$ 

Using Eq. (|73j) . Eq.lfTHjl. and the result the contribution of Eq. (JHH|) to E$ for p may 

be easily computed 

n M x 



< M\H q ^ M - H q \M >< M\H, 



3N r - 13 N, 



H q \M > 



g 6 ad 



c 



N 



c 



q ( tM 

M — e/ ' p 

N 2 -18iV c + 37 N c ~ 
(N c -3) 2 2 



(104) 



Collecting all these contributions, one finds that the iV-dependent terms cancel against 
the quark Hamiltonian contribution to Eq^ [see Eq.([47|)] and the remaining terms are N- 



25 



independent. This is a very good check of our computation since the meson mass should be 
an intensive quantity. 

Adding up all the iV-independent terms the contribution due to the quark Hamiltonian 
is given by 

(4) _ „(4) 31180 - 48674iV c + 29051iV£ - 8183iVg, + 1069iV£ - 5LZ\% 

2g«aCl{N c -±){N c -mNc-2) ■' ,) 



where is the contribution to the fourth order vacuum energy due to the sole quark 
Hamiltonian l|47p. The results for the other mesons are listed below (see also Appendix D 
for the results for each matrix element): 

(106) 
(107) 
(108) 
(109) 
(110) 
(111) 

Equation (jBT)) yields the magnetic contribution to E$ and it is at the second order in the 
strong coupling expansion. This term gives the same result for all the mesons and reads 

(E^Uagnetic = + g 6 a ( 2 N* - 1)(N 2 - 1) + g & a(2N c - 3)(iVj - 1) ^ 

Again the iV-dependent term in Eq. ()1120 vanishes in the difference between Eq. (11120 and 
the magnetic contribution to the fourth order correction to the vacuum energy E^, Eq.(|470. 
As for the second order, the fourth order correction to the meson energy is finite and well 
defined after rescaling the coupling constant g 2 Nc —>■ g 2 and taking the 't Hooft limit (large 
N c with g 2 N c fixed). 



4 4) 


P (4) 


29452 - 


45650iVc + 29963JVg - 7471iVg + 949A^ - 43iV£ 




2 9 eaC$(N c - 4)(N C - 3) 3 (iV c - 2) 




P (4) 


29452 - 


45650iV c + 26963A^ - 7471A^ + 949iV£ - 43iV£ 




2g^aCl(N c - A)(N C - 3) 3 (iV c - 2) 


Eif 


P (4) 


36172 - 


56386iV c + 33803JV£ - 9671iV"£ + 1309iV£ - 67iV"£ 




2g*aC${N c - 4)(iV c - 3) 3 (iV c - 2) 


E^ ] 

ai 


P (4) 
— ^q0 — 


124972 - 


- 205002iV c + 130895JVg - 40619JV& + 6121JV£ - 359A^ 




2g°aCi(N c - 4)(iV c - 3) 3 (iV c - 2) 


P (4) 
E h 


F (4) 


29452 - 


4565(W C + 26963A^ - 7471A^ + 949JV£ - 43iV£ 




2g e aCl{N c - 4)(N C - 3) 3 (N C - 2) 




_ P( 4 ) - 
^q0 — 


76396 - 


127034AT C + 8311 liV£ - 26803iV£ + 4273iV£ - 271iV£ 




2g*aCi(N c - 4)(N C - 3)%N C - 2) 
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D. The meson spectrum 



The lattice excitation masses are given by subtracting the energy of the ground state 
from the energies of the excitations 

m M = E M — E 

Using the results of Eqs. (J52j) - (J55j) and of Eqs. (jlU5j) - (jll2j) . one gets in the large N c limit 

m W0 = ^(i + 6e 2 -171e 4 ) (113) 
a 2 1 

m p = ^-(- + 6e 2 - 203e 4 ) (114) 
a 2 1 

m u = — ( 7 + 6e 2 - 171e 4 ) (115) 
a 4 

a 2 1 

m bl = — (- + 10e 2 - 267e 4 ) (116) 
a 4 

a 2 1 

m a = — ( - + 14e 2 - 1435e 4 ) (117) 
a 4 

m f = ^_(I + i4e 2 - 875e 4 ) (118) 
a 4 

a 2 1 

m/o = £_(_ + I8e 2 - 1083e 4 ) (119) 

where e = 1/g 2 and with (7 we indicate the rescaled coupling constant g 2 Nc —>■ g 2 ■ Equa- 
tions (jll3|) - (jll9|) provide the value of the meson masses up to the fourth order in the strong 
coupling expansion. 



V. LATTICE VS. CONTINUUM 

The series given in the preceding section are derived for large g 2 . Since from renormaliza- 
tion group arguments for an asymptotically free theory g 2 = —c/ In a for small a, the series 
for the meson masses are valid only for large lattice spacings. To compare the results of the 
strong coupling expansion with the continuum theory one needs some method of continuing 
the series to the region in which g 2 = 0, i.e., e = 00. To make this extrapolation possible, 
it is customary to make use of Pade approximants [2^, which allows one to extrapolate a 
series expansion beyond the convergence radius. For this purpose one should consider the 
mass ratios, expand them as power series in y = e 2 and then use [1, 1] Pade approximants 
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by writing the mass ratios in the form 

pi = 1 + a V 
1 1 + by 

where a and b are determined by expanding to order y 2 and equating coefficients. In the 
continuum limit this ratios yields a/b. 

The results obtained with this method for the mass ratios of the single link mesons are 
listed in Table 1. For each mass ratio considered here, the [1, 1] Pade approximant exists for 

TABLE I: Mass ratios of the single-link mesons 



Mass ratios Our results Banks et al. Experimental values 



my 
m bl 


0.75 


0.86 


0.11 


m bl 


0.75 


0.86 


0.63 


m p 
m bl 


0.71 


0.81 


0.62 


m bl 
m fo 


0.82 


1.03 


0.90 


m bl 
m ai 


0.95 


0.93 


0.98 


m bl 
m f 2 


0.92 


1.00 


0.97 



positive values for a and b. Therefore, the extrapolation from y = to y = oo is singularity 
free in our approximation. 

We considered the ratios between meson masses involving the b\ meson. For these ratios 
the results are in very good agreement with the well known experimental values except for 
the pion mass. This foreseeable failure is due to the lack of full chiral symmetry in the 
theory for large lattice spacing. The n-p splitting is so tiny because of the lack of significant 
spin-spin forces in the first four orders of strong-coupling perturbation theory. Magnetic 
field effects, loops of flux, are just not important through this order. However, at sixth and 
eighth order such effects are important and should provide improved results. 

The mass of f seems to be quite large but in agreement with the experimental value of 
the meson / (1370). 
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Our results show that the large Nc expansion provides a very good and systematic 
theoretical setting to evaluate physical quantities of phenomenological interest. 

VI. CONCLUDING REMARKS 

In this paper we used the strong-coupling expansion of lattice QCD to compute -for 
large Nc- the low-lying unflavored meson spectrum. Our large Nc Hamiltonian approach 
with staggered fermions evidences that the possible ground states of strongly coupled lattice 
QCD are those of a spin Ac/2 antiferro magnetic Ising model. Choosing one of the two 
ground states amounts then to the spontaneous breaking of the discrete chiral symmetry 
corresponding to translations by a lattice site. As a consequence a nonvanishing chiral 
condensate should arise; it would then be interesting to compute -within our formalism- the 
chiral condensate on the lattice and then compare it with the results of numerical simulations. 

Mesons are created by operators that, acting on the vacuum, create qq states with the 
desired quantum numbers. Their energy is computed up to the fourth order in the strong 
coupling expansion. The meson masses are obtained by subtracting the vacuum state energy 
from the energies of the excitations. After rescaling the coupling constant g 2 Nc — > g 2 , 
according to the 't Hooft prescription, and taking the Nc — > oo limit one finds that the 
series for the meson masses are well defined. Since it is expected that the continuum limit 
occurs without any phase transitions, the series are analytically continued by using Pade 
approximants. 

With the exception of the pion, which turns out to be degenerate in mass with u, the 
results we obtained are, already at the fourth order, in good agreement with observed values. 
Higher orders of the strong coupling expansion should evidence the splitting between the 
masses of tiq and u since the effects of the magnetic field would be more pronounced. 

The strong-coupling limit of QCD is not universal, thus adding an irrelevant operator to 
the Hamiltonian leads to the same physical predictions in the continuum limit. This allows 
one to introduce arbitrary parameters, the coefficients of these irrelevant operators, which 
are then fixed by fitting the experimental data. We have shown that, in the 't Hooft limit, the 
meson masses can be computed without introducing irrelevant operators and consequently 
without arbitrary parameters. Irrelevant operators are usually introduced Q], in fact, 
mainly because in this way the meson masses are well defined for Nc = 3. Our large Nc 



29 



approach does not need any irrelevant operator; it yields series expansions for the meson 
masses which are free of divergences and thus well defined. 

We determined the ratios between two single link mesons even if these ratios are expected 
to be less reliable than the ratios involving meson masses relative to the nucleon mass. The 
reason is that all the mesons are degenerate at zeroth order so small differences of second 
and fourth Taylor series coefficients control the [1, 1] Pade approximants, and a considerable 
amount of information in the series is lost. The good agreement between our results and 
the experimental values for the ratios between meson masses shows that the large Nc limit 
is very effective in the strong-coupling region also. It should be interesting to extend our 
procedure to provide a large N c evaluation of the ratios between the meson masses and the 
nucleon mass. 
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APPENDIX A: THE STAGGERED FERMIONS FORMALISM 



In this appendix the well known formalism for staggered fermions originally 



by Kogut and Susskind 



2l\ is reviewed. Other formulations exist in literature 



13 



developed 



221 but in 



this paper we use the notation and the conventions of refs. 

In this approach the two nonstrange quark fields u and d are represented by a single- 
component lattice field ^(r). In a three-dimensional cubic lattice with sites labeled by 
triplets of indices r = (x,y,z) one may define the lattice derivative operator 



Vj^(r) = — ip{f + i) — ip{f — i) 
With the definition (jAl)) the massless Dirac equation 



(Al) 



if) = —a ■ Vijj 



becomes, on the lattice, 



tp{r) = V] a ■ h[ip(r + h) — ip(r — n)] 



(A2) 
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FIG. 1: Labeling of lattice sites 
In Eq. ()A2|) a is the vector of Dirac matrices in the representation 





d \ 


a = 


fo 




7o = 






:) 




\° - 1 ) 









(A3) 



The dispersion law for Eq. (|A2|) is given by 

uj = (ce x sin l x + a y sin l y + a z sin i a ) / a 

and the low-energy spectrum is found at the eight corners of a cube in k space. 

It is customary to remove the degeneracy by reducing the degrees of freedom. One may 
consider a one-component field, <f>(r), per lattice site which will represent a component of 
the Dirac filed. The lattice is subdivided into four sublattices to place the four components 
of a conventional Dirac field and the corners of a unit cube are labeled as shown in Fig.l. 
By further subdividing the lattice in sublattices there are two complete and independent 
Dirac fields which exhaust the low-frequency spectrum of the Dirac equation. The lattice is 
now divided into a sublattice of sites for which y is even and another of sites for which y is 
odd. The fields are conveniently relabeled as follows 

A = ft (y = even) 
A = 9i (y = odd). 

Writing Eq. (jA2|) in Fourier transformed variables, one finds for / and g 

uiaf = (a z sin k z + a x sin k x )f + (a y sin k y )g (A4) 
uag = (a z sin k z + a x sin k x )g + (a y sin k y )f. (A5) 
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There are two combinations of / and g which are conventional Dirac fields for long wave- 
lengths 

u i = f l + g l (2 = 1,2,3,4) (A6) 

and 

di = h ~ 92 
d 2 = -(/i - 9i) 
d 3 = ~(fa ~ 9i) 

di = h- 93 (A7) 

This is easily seen by summing up Eqs. flA4|) and (|A5|) . One finds for the component of the 
u field in the representation (jA3|) 



auuui(k) = sin k x u± — i sin k y u^ + sin k z u 3 
auiu 2 (k) = sin k x u 3 + % sin k y u 3 — sin k z u^ 
au)Uz(k) = sin k x u<i — i sin k y ui + sin k z u\ 

auiUi(k) = sin k x U\ + % sin k y U\ — sin k z u 2 (A8) 

— * 

and identical equations for the ef s. If one now considers those normal modes with k small 
and set k = aK, Eq. (|A8|) . for small k, takes the form 

ujui(K) = K x u± — iKyUi + K z u 3 
uu 2 (K) = K x u 3 + iK y u 3 - K z Ui 
uju 3 (K) = K x u 2 - iK y u 2 + K z u\ 

uju a (K) = K x u\ + iKyU\ — K z u 2 (A9) 

and the same with u < — > d. Equations (jA9|) are the normal Dirac equations for u and d in 
the representation (|A3|) . 

One may now turn to the free massless Dirac Hamiltonian in the continuum 

H = —i(v) { OLidiU + SoLidid) 

and write it in the lattice form by means of the definitions (|A6|) and (jA7|) for the u and d 
fields in the representation (jA3|) . In terms of the one component field 0(f) one gets 

H = -f E [fa f (*W+ - h. c .){-if + y 



2a 

r 
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+ ($ (r)(j)(r + x) — h.c.) 

-%{(j)\r)(j){r + y) + h.c.){-l) x+y ] . (A10) 

Equation (|A10J) may be written in a more symmetric form by defining the following 
functions on the lattice 

D(x, y) = l -[{-lf + {-If + + 1] (All) 

A{ n ) = -^=\i n - 1 ' 2 + (-l) n ^ 2 ]. (A12) 
v2 

On the lattice sites D and A are either 1 or -1 and therefore 

D 2 = A 2 = 1. 

Furthermore, D and A satisfy the relations 

D(y,x)D(y,x + l) = (-iy 

D(y + l,x + l)D(y,x) = (-iy + y +1 

A(y)A(y + l) = (-iy. (A13) 

Defining the field 

*(f) = (-l)»(i) s+ *A(3/)Z>(a:,z)0(r) (A14) 
and using the previous definitions for A and D the lattice Hamiltonian becomes 

H * = \Y, [^(r + x)^(f)(~iy 

r 

+ y)V(r)(-l) x 
* f (f + i)^(f)(-l) y + /i.e.] , (A15) 

which is the quark Hamiltonian used in our paper. 



APPENDIX B: CONSTRUCTION OF AN EIGENSTATE OF H e 

In this Appendix we explicit the construction of eigenstates of the unperturbed Hamilto- 
nian H e in order to evaluate the function of H e in the perturbative expansion. Using Eq.(JHJ), 
Eq.© and Eq.(j22J), one finds 

[H e ,U{ B [r,n]Ul D [r>,m}]\0> = £ {C 2 U\ B [r, h]rf OD [P, m] 
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^^U\ B [f,n}Ul D [f,n]5([r,h] - [r',m])) |0 > .(Bl) 



Using Eq. ()Bl|) . one finds 



[H ej UUr,n]Ul D [f,h]}\0 >= g* {Na 2) U\ B [r, h]tf CD [r, n]jO > . (B2) 

One may now look for an eigenstate of H e with eigenvalue g 2 C2/a 
H e (u\ B [T,h]Uh D [r\m] + aUi B [r,n]tf CD [r,n}) |0 



>= 

^C 2 (u{ B [f, h]Ul D [J, m] + aU\ B [r, h]tf CD [r, h]) |0 > (B3) 

Taking into account Eq. (|B1J) and Eq.(jB2j) one gets a = — 1 and then the pertinent eigenstate 
is 



He (u{ B [r,n}Ul D [ri,m] - U{ B [f, npl D [f, n]) |0 >= 

(^ B [r, n]4 D [rt, m] - U\ B [f, n]C/J D [f, n]) |0 > (B4) 

Using Eq. ()B4|) . one may evaluate the function of H e appearing in the perturbative expansion 

n 



f(H e 



Ei 0) - H e 



From Eq. (|B2|) . one has 



f(H e pUr,np j CD [r,n}\0 >= - -U{ B [r, np^ CD [r, n]|0 > . (B5) 

g {i\c + iJiiVc — 

Using Eq. (|B5|) . one gets from Eq. (jB4|) 



/(iy e )C/I B [r,n]^ D [r',m]|0 >= -— (^[^feK™] 

^U\ B [r,hpl D [r,h]5([f,h] - [P,m])) |0 > ■ (B6) 



APPENDIX C: INTEGRATION OVER SU(N C ) 

In this appendix a table of the integrals over the group elements of SU(Nq) needed in 
the paper is provided. 

It is well known that a basic ingredient to formulate QCD on a lattice is to define 
the measure of integration over the gauge degrees of freedom. Unlike the continuum gauge 
fields, the lattice gauge fields are SU(N C ) matrices with elements bounded in the range [0, 1]; 
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Wilson [tJ proposed an invariant group measure, the Haar measure, for the integration over 
the group elements. The integral is defined so that, for any elements g% and g<i of the group, 
one has 

J dUf(U) = J dUf(U 9l ) = J dUf(g 2 U), (CI) 

with f(U) a generic function over the group. When used in nonperturbative studies of gauge 
theory, the definition (jCl|) avoids the problem of introducing a gauge fixing, since the field 
variables are compact. The measure is normalized as 

dU = 1. (C2) 

The strong coupling expansion for an SU (Nc) gauge theory depends on the following iden- 
tities for integration over link matrices 



dUU ab = J dUU ] ab = (C3) 
dUU ah U\ d = -UoAc (C4) 

dUU ab Ul d U ef Ul h = 2 _ (SadSbcSehSfg + $ah$bg$cf$de) 

c 

1 

- (5ad5bg5eh5fc + 6 ahhJ> ed& '} g) ■ (C5) 



N C (N 2 C - 1) 

One also needs the group integral over six elements, which occurs at the fourth order in the 
strong coupling expansion of the mass spectrum. The pertinent integral is 

r : N 2 — 2 

/ dUUahU^U e fU h UiiU mn = 2 — c 2 -rz(5 a dS e h5i n 5b C 5fg5i m + S a dS en 5hiSb C 5f m 5 g i 

J iVcVJV c — ij(iv c — 4J 

+^ah^de^inhg^cf^lm + S an 5diS e hSbmS c l5 f g + 5 an 5de^hi^bm^c}^gl + & ah$ di& en$bg& cl& fm) 

l^J\f 2 ~£j(dadfien8hi8bc$ fg$lm + ^ah^de^in^bc^ fg^lm + ^an^di^eh^bc^ f g$lm 

+3ad3ehd in^bc^fm^gl + ^an^de^hi^bcO fmOgl + 5ah&di&enSbc& fm^gl + ^ad^eh^in^bg^cfSlm 
+5an5de5hi5bg5 c f5lm + ^ah^di^en^bg^cf^lm + S a dS e hSi n SbmS c l^ fg + S an 5d e ShiSbmS c lS fg 
+OahOdiO e nUbmV c lOfg + S a dS en 5hiSbm^cfSgl + ^ahOdeOi n Obm^cf^gl + ^an^di^ehObm^cf^gl 
+^ad^en^hi^bg^cl^ fm + ^ah^de^inhg^cl^ fm + ^an^di^eh^bg^cl^ fm) 

+5an5di5 e h8bc8fm8gl + S a dS en 5hiSbgS c f5l m + S an SdiS e hSbgS c fSl m + SadbenShiSbmbclbfg 
+8ah8de8in8bm8cl8fg + ^ad^eh^in^bm^cf^gl + ^ah^di^en^bm^cf^gl + S a dS e hSi n 5bgS c l5 f m 
+SanSdeShiSbgSdSfm)- (C6) 
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APPENDIX D: MATRIX ELEMENTS 



In this appendix the matrix elements useful for the computation of the n, u, b\, a±,f2 
and /o energies at the fourth order in the strong-coupling expansion are reported. 

I. 7T 



< M^\Hlk M Hlk M H q k M H q \M >= 
18(7V C -1) 4(iV c -l) 2 



(N c - 1) 
N c -2 



-N C N - 58 
2 ° 



4(N C - 1)- 
" (N c - 3) 2 



iVc-3 



(iV c iV - 22) 
iVc-2 



+ 



6 



N c -A 



N C /2N + 2N C - 14 
iVc-2 



+ 



6 



7V c -4 



(Dl) 



< M\HlA M HlA M H q A M H q \M^ > 
, 12(7V C - l) 2 



iVc-2 



+ AM 87V C + 



47W h72 - 12iVc 
TV 2 7V c (iV c -l) 2 \ _ 

2 ~ 2(iV c -2) ) 4 



/V 2 



(D2) 



< M\H q K M Hlk M H q K M Hl\M^ > 



+N Un c + ^ 



(N c - l) 4 



(N c -2)(N c -3y 



-33 - 127V C - 



12 + ^7V 



N 2 

—!=-N 2 



(D3) 



< M\Hlk M H q k M Hlk M H q \M^ >-- 

12(N C -1) 2 , T ( AAT N c 

H — — + N 4N C H - 

iVc-2 ^ 6 2 



4^ 6 aC| 1 
iV c (iV c -l) 2 ' 

2(iV c -2) , 



12AT C 

N 2 
-^C N 2 



(D4) 



< M\H q k M Hlk M H\k M H q \M^ >= 

*t(**t N c N c (N c -lf 
+N 4N C + — 



Ag^aCl 



-32 - 12iV c + 



< M\H q k M Hlk M Hlk M Hl\M > 



2 2(N C - 3)(N C - 2) J 8 
1 



AT 2 

Nc N 2 



4g*aCl 



14 + AN C - N C N - 



12(N C - l) 3 
(N c -3)(N c -2) 

(D5) 

4(iVc - l) 3 
(N c -3)(N c -2) 
(D6) 



< M\H q k M Hlk M H\k M Hl\M >= 



12 + AN C - 4( ^ _ - N C N 



(D7) 
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2. uj 



n 



M 



-H q \M >< M\H q 



n 



N c - 18N C + 37 



9 6 CI 



'M 

3N C 



N c -3 2 



(n c - sf 



H q \M > 



N c 



(D8) 



< M^\H\k M Hlk M H q k M H q \M >= 

-J— J -96 + ("-Nc) N-^N 2 - ffg " 1} 
^aCf \ + U C 7 4 iV c - 2 

18(N C - 1) 4(iV c - l) 2 \(N C N - 22) 



-NqN - 58 



iV c -4 
4(iVc - l) 2 
' (N c - 3) 2 



N c -3 



N 



+ 



6 



N 



c 



N C /2N + 2N C - 14 
iVc-2 



+ 



6 



N c -A 



(D9) 



< M\HlK M Hlk M H q K M H q \M^ >-. 



1 



4£ 6 aCf 



-72 - 12iV c + 



12(N C - if 



Nr 



,r ,r N C N C (N C -1) 2 \ N C ~ 

+N 8N C + — r~ -N 2 

V 2 2iV c -2 J 4 



(DID) 



< M\H q A M HlA M H q A M Hl\M^ >= 



N C-N 2 



(N c - l) 4 



N 2 C 



-41 — 12Nc + N \ 4Nc + 



8 



(N c -2)(N c -3f 



-12 + ^yN 



(Dll) 



< M\HlK M H q K M H\K M H q \M^ >-- 
° 2 2(iV c - 2) , 



4<? 6 aC 2 3 



-42 - 12iV c + 



12(N C - if 
N c -2 



(D12) 



< M\H q k M Hlk M H\k M H q \M^ >= 
( N 2 



iV 2 7V c (iV c - l) 3 \ 



-36 - 12N C + 



< M\H q K M H\K M HlK M Hl\M > 



2 2(N c -3)(Nc-2)J 8 
1 



N °N 2 



22 + 4N C - N C N 



12(N C - l) 3 
(N c -3)(N c -2) 

(D13) 

4(iVc - l) 3 
(N c -3)(N c -2) 
(D14) 



< M\H q k M Hlk M Hlk M Hl\M >-- 



24 + AN C - — — - NcN 



N c -2 



(D15) 



37 



3. 6i 



n 



< M\H q ^ M H q \M >< M\H q 

<J M ~ tie 

-3N C - 13 



n 



M 



Hg\M > 



N c ' 
-N 

N c -3 2 



N c - 18N C + 37 , Nc 



(N c - 3) 2 



+ -N 



(D16) 



< M^\Hlk M Hlk M H q k M H q \M >= 

-J— { -76 + f^iV^ iV - Utrf - ^ " 1} 
Ag^aCl\ V2 V 4 iV c - 2 

18(iV c - 1) 4(iV c - l) 2 \(N C N - 22) 



-N C N - 58 



iV c -4 
4(iVc - l) 2 
' (N c - 3) 2 



N c -3 



N 



+ 



6 



c 



N 



c 
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